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Evidence for dark matter self-interactions has recently been reported based on the observation of 
a spatial offset between the dark matter halo and the stars in a galaxy in the cluster Abell 3827. 
Interpreting the offset as due to dark matter self-interactions leads to a cross section measurement of 
cdm/ui ~ (1 — 1.5) cm^ g~^, where m is the mass of the dark matter particle. We use this observation 
to constrain singlet scalar dark matter coupled to the Standard Model and to two-Higgs-doublet 
models. We show that the most natural scenario in this class of models is very light dark matter, 
below about 0.1 GeV, whose relic abundance is set by freeze-in, i.e., by slow production of dark 
matter in the early universe via extremely tiny interactions with the Higgs boson, never reaching 
thermal equilibrium. We also show that the dark matter abundance can be established through 
the usual thermal freeze-out mechanism in the singlet scalar extension of the Yukawa-aligned two- 
Higgs-doublet model, but that it requires rather severe fine tuning of the singlet scalar mass. 


I. INTRODUCTION 

Understanding the nature of dark matter is one of the 
holy grails of modern particle physics. Up to now, dark 
matter has been observed only through its gravitational 
interactions with Standard Model (SM) particles (for a 
recent review, see, e.g.. Ref. [T]). 

Very recently, the first evidence for dark matter self¬ 
interactions has been reported [2] based on observations 
of four elliptical galaxies in the inner 10 kpc core of galaxy 
cluster Abell 3827. Using gravitational leasing. Ref. [2] 
reconstructs the dark matter halos of the four galaxies 
and observes that at least one of the halos is spatially 
offset from its stars by a distance of A = 1.62to;49 kpc. 
Dark matter self-interactions can lead to such an offset: 
the dark matter halos of the infalling galaxies will be 
impeded as they pass through the dark matter associated 
with the cluster. The magnitude of this effect scales with 
the self-interaction cross section, which can therefore be 
extracted from the size of the offset. 

Kahlhoefer et al. [3] consider the case of dark matter 
contact interactions, which will be relevant for the scalar 
dark matter models we consider below. Whereas the ef¬ 
fect on the halo of frequent dark matter interactions can 
be described as a drag force, in the case of rare, high- 
momentum-transfer contact interactions, this picture is 
inaccurate. Rather, in the contact-interaction scenario 
the observed separation is due to the evaporation of dark 
matter from the halo. Most dark matter particles in the 
halo will not interact, but those that do scatter will be 
ejected from the halo in the backward direction. Even¬ 
tually the scattered particles will become substantially 
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separated from the halo, but at short time scales, they 
remain observationally associated with the galaxy. This 
leads to a deformation of the halo and an offset between 
the centroids of the galaxy and the halo. Determination 
of the size of the offset is fairly sensitive to how the cen¬ 
troid is defined, but Ref. [5] finds that a cross section of 

cdm/to ^ 1-5 cm^ g“^ (1) 

could account for the offset observed in Ref. [2]. This 
is in slight tension with results from previous studies of 
other astrophysical systems EHU, which find an upper 
bound on the cross section of 

ctdm/to < 1 cm^ g“b (2) 

We will use Eqs. ([^ and ([^ to set a range for the esti¬ 
mated self-interaction cross section. In particle physics 
units, this corresponds to 

ctdm/to - (4.7 - 7.0) X 10^ GeV^ (3) 

Of course, it is possible that the offset may be an as¬ 
trophysical artifact. Nevertheless, if a dark matter self¬ 
interaction interpretation holds up, it would represent a 
profound advance in our knowledge of the nature of dark 
matter. It is therefore of great importance to explore its 
implications for models of dark matter. 

In this paper we consider models in which the dark 
matter is a gauge-singlet real scalar particle S that inter¬ 
acts with the rest of the SM through the “Higgs por¬ 
tal” via an operator of the form where $ is 

the SU(2)i-doublet SM Higgs field. The singlet scalar 
is kept stable by imposing a symmetry S —)■ —S on the 
Lagrangian. In Sec. |ll] we consider the minimal model 
in which a single gauge-singlet real scalar is added to 
the SM [5H25]. and show that it cannot explain the ob¬ 
served dark matter self-interaction cross section once we 
impose the requirement that the annihilation cross sec¬ 
tion of the singlet scalars into SM particles yields the 
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correct relic abundance from thermal freeze-out together 
with the constraints on the invisible branching fraction of 
the Higgs boson from the CERN Large Hadron Collider 
(LHC) [2H1 ES] ■ The correct relic abundance can, how¬ 
ever, be achieved in this model through thermal “freeze- 
in,” in which the dark matter is produced slowly in the 
early universe through extremely weak interactions with 
the SM Higgs boson [T31[inil3I]- In Sec. HI we attempt to 
salvage the usual thermal freeze-out picture by consider¬ 
ing the real singlet scalar extension of two-Higgs-doublet 
models (2HDMs) [32H32]- We find a small viable region 
of parameter space in which the singlet scalar mass is 
tuned to be very close to half the mass of the second, as- 
yet-undiscovered CP-even neutral Higgs boson of these 
models. In Sec. m we conclude. 


II. MINIMAL SINGLET EXTENSION OF THE 
STANDARD MODEL 
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We write the scalar potential for the SM plus a real 
singlet as 

V = -h ^SS + -b ApS'S'$I$ -b AsS”^, 

(4) 

where we have imposed invariance under S' —S' to keep 
the singlet stable. The physical SM Higgs mass nift ~ 
125 GeV [33] and the SM Higgs vacuum expectation value 
(vev) V = ~ 246 GeV |33| fix two of the 

scalar potential parameters, 

= XhV^ ~ (88.4 GeV)2. (5) 

For the potential to be bounded from below we require 
that As > 0 and \p > —2y'XhXs- To keep the singlet 
from getting a vev and thereby ceasing to be stable, we 
require that 

(6) 

where ms is the physical mass of the singlet scalar in the 
electroweak-breaking vacuum. This implies that at most 
one of /x| or Ap can be negative. If ^| < 0 the potential 
can develop a second minimum in which the singlet gets 
a vev and the doublet does not; to avoid this possibility 
we require that > 0. 

The Feynman rules involving the singlet that will be 
important here are 

SSSS : -24iXs, 

SSh : —2iXpV. (7) 

We will consider dark matter scattering SS SS in 
three mass ranges. We work in the minimal real singlet 
scalar extension of the SM, but what follows will inform 
the calculations in the extended Higgs sector that we con¬ 
sider later. The relevant Feynman diagrams are shown 
in Fig. [2 


FIG. 1. Feynman diagrams for the SS — >■ SS scattering pro¬ 
cess in the real scalar singlet extension of the SM. 

A. Heavy dark matter, ms ^ mu 

When S is very heavy compared to the SM Higgs mass, 
the Higgs-exchange diagrams in Fig. are suppressed by 
momentum factors compared to the four-point diagram, 
so we neglect them for simplicity. The scattering cross 
section per unit dark matter mass can then be written as 

, 9A| 

Perturbative unitarity of the 2 —> 2 scattering ampli¬ 
tude for SS ^ SS at high energy sets an upper bound on 
Ag. Requiring that the zeroth partial wave amplitude Og 
satisfy |Re ao| < 1/2 and taking into account the identi¬ 
cal particle factors of 1/ \f2 for the initial and final states 
yields 

As < 27r/3. (9) 

Unfortunately, after imposing this constraint we find that 
the heavy-singlet scenario is inconsistent with the mea¬ 
sured cross section in Eq. Indeed, reproducing the 

observed dark matter self-interaction cross section would 
require ms ^0.1 GeV; we will consider this low-mass 
scenario below. 


B. Dark matter near the Higgs resonance, 

ms ~ mh/2 

When the scattering is near resonance, the s-channel 
Higgs exchange diagram in Fig. dominates. Defining 
the hSS coupling 

gssh = -‘2XpV, 


( 10 ) 
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we can write the hSS Feynman rule as igssh- The s- 
channel Higgs exchange diagram in Fig. then yields a 
dark matter self-interaction cross section per unit mass 

00 

o’dm/ws 1287rm|.(4TO| — 


IdTrmKms — mh/^Y ’ 

where in the last expression we have taken ms — mu!^- 
Perturbative unitarity of the four SS —> scatter¬ 

ing amplitudes, where ipi are the four real scalar compo¬ 
nents of the SM Higgs doublet constrains 



FIG. 2. Feynman diagrams for the annihilation of singlet 
scalars into SM particles in the real singlet scalar extension 
of the SM. 


|Ap| < dTT, or |gss/,| < St:v. (12) 

Achieving the observed dark matter self-interaction cross 
section requires that ms be severely tuned such that 
|ws — m/j/2| < 0.11 GeV, where we use the lower value 
of ctdm/tos in Eq. © to obtain an upper bound. This 
represents a tuning of about two parts per mille. 

We now show that this scenario cannot be reconciled 
with the observed dark matter relic density because the 
annihilation cross section into SM particles through the 
Higgs resonance is several orders of magnitude too large. 


1. Relic density 


The annihilation cross section required to obtain the 
correct thermal dark matter relic density for ms ~ nr;j/2 

is [15] 

av^e\ ^ 2.2 X 10-2® cm® s-^ ~ 1.9 x 10"® GeV^. (13) 


For ms ~ myj^^ the singlet scalars annihilate into 
SM particles predominantly through s-channel Higgs ex¬ 
change, as shown in Fig. The temperature of the 
thermal bath at the time of freeze-out for a dark matter 
particle with the correct relic density is T ~ msl‘2.Q ^ 
5 GeV [T]. This is much larger than the mass split¬ 
ting \ms — mh/‘2.\ < 0.11 GeV required to obtain a large 
enough dark matter self-interaction cross section, so ther¬ 
mal effects during freeze-out cannot be neglected. 

The annihilation cross section can be written as 


CV^el 


9„2 p^iSM 
■^ySSh^ tot 


(14) 


where Ecm = ‘2ms + insv^ei/^ is the collision center-of- 
mass energy, Vrei is the relative velocity of the two dark 
matter particles in the center-of-momentum frame and 
~ 4.1 MeV is the SM Higgs total width [46] . 


^ The contribution to the center-of-mass collision energy of the 
present-day kinetic energy will be totally negligible compared to 
the mass splitting that we will find, so we neglect it here. 


The singlet scalar can be heavier or lighter than half 
the Higgs mass. For a singlet scalar heavier than half 
the Higgs mass, the thermal motion pushes the anni¬ 
hilation process farther from resonance. Imposing the 
self-interaction cross section measurement and integrat¬ 
ing the annihilation cross section over the Boltzmann dis¬ 
tribution with T ~ tos/20 yields 

(o-^^rei) > 0.03 GeV-2, (15) 


where the bound is saturated for |gss/i| as large as pos¬ 
sible. This is many orders of magnitude larger than 
the required cross section for thermal freeze-out given 
in Eq. (13), resulting in a singlet scalar relic abundance 
way too small to account for the observed dark matter. If 
the singlet scalar is lighter than half the Higgs mass, the 
thermal motion pushes the annihilation process onto the 
resonance resulting in an even larger thermally-averaged 
cross section and even smaller relic abundance. These 
large cross sections are a consequence of the large SSh 
coupling and the tuning of the singlet mass close to the 
Higgs resonance that is required to achieve the observed 
dark matter self-interaction cross section. 


C. Low-mass dark matter, ms <C mh 


We now consider the situation of low-mass dark mat¬ 
ter. In this case all four diagrams in Fig. [^ contribute. 
At low center-of-mass energies, the Higgs boson h can be 
integrated out, yielding an effective four-S' coupling with 
Feynman rule —24iAeff, where 


Aeff = As — 


9ssh 
8ml ■ 


(16) 


SS ^ SS scattering in the low-energy effective theory 
must also satisfy perturbative unitarity, so as in Eq. © 
we obtain 


|Aoff| < 27r/3. 


(17) 


The dark matter scattering cross section per unit mass 
as in Eq. (|^ becomes 




( 18 ) 
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FIG. 3. Favored region in ms and the effective 45" conpling 
Acff defined in Eq. ( |16| l for a dark matter self-interaction cross 
section in the range 1.0-1.5 cm^ g“^. Values of |Aeft| above 
27r/3 ~ 2.094 (horizontal dashed line) are excluded by per¬ 
turbative unitarity considerations. The vertical dotted lines 
show the thresholds for SS annihilation into pairs of electrons, 
muons, and neutral pious. The charged pion threshold is at 
the right edge of the plot. 


The range of ms and AeS consistent with the cross section 
measurement from Ref. El [Eq. @] is shown in Fig. 

In particular, the singlet scalar dark matter must be very 
light, with mass below about 0.1 GeV. 


1. Invisible Higgs decays 


One can also derive an upper bound on gssh from the 
measured signal strengths /ti of the Higgs boson in var¬ 
ious SM decay channels. In the singlet scalar extension 
of the SM, these signal strengths are all modified by the 
common factor 


ft* = M = 


p/iSM 

^ tot 




V{h^ SS) 


( 22 ) 


Combining all channels, ATLAS [17] and CMS [15] re¬ 
spectively measure an overall signal strength of 

= (ATLAS), 

^ = (CMS). (23) 

In the real singlet extension of the SM, /r < 1, therefore, 
relative to this best-fit point, the = 4 lower bounds 
on the signal strength are 

fi > 0.85 (ATLAS), 

fi>0.74 (CMS). (24) 

These in turn lead to upper bounds on gssh^ ot’ equiva¬ 
lently |Ap|, of 

glsh < 9-1 GeV^ or |Ap| < 6.1 x lO'^ (ATLAS), 

< 18 GeV^ or |Ap| < 8.6 x lO'^ (CMS). (25) 


2. Relic density: thermal freeze-out 

The annihilation cross section required to obtain the 
correct thermal dark matter relic density varies slightly 
with ms [45] . For ms around 0.1 GeV, the annihilation 
cross section needs to be [IS] 


When ms ^ my, the non-observation of invisible 
Higgs boson decays constrains Ap. 

The partial width for h ^ SS is 


r{h SS) 


9ssh 

327rTO/i 



9ssh 

327rTO/i ’ 


(19) 


where the last expression holds for ms ^ mh- The in¬ 
visible Higgs branching fraction is defined as 


BR(/i —)■ invis.) 


T{h SS) 

FfsM+r(/i^55)' 


( 20 ) 


The strongest constraint on BR(/i —invis.) currently 
comes from the ATLAS search for invisible Higgs decays 
in vector boson fusion production, with a 95% confi¬ 
dence level upper limit of BR(/i —>■ invis.) < 0.29 for 
a 125 GeV Higgs boson, assuming a SM production cross 
section [29]. This leads to an upper bound on gssh: 
equivalently |Ap|, of 

glsh < 21 GeV^ or |Ap| < 9.3 x 10"^ (21) 


This constraint is very strong compared to the perturba¬ 
tive unitarity constraint |Ap| < 47r. 


crurei ^ 4.8 X 10"^® cm® s"® ~ 4.1 X 10"® GeV"®. (26) 


In the low-mass range, the singlet scalars annihilate 
into SM particles predominantly through s-channel Higgs 
exchange, as shown in Fig. [^ This annihilation cross 
section can be conveniently expressed for ms 'C m^ as 


CTUrel 


mf ms 


(27) 


where = 2ms) is the would-be total width of 

the SM Higgs boson if its mass were equal to 2ms- The 
largest possible values of av^ei are obtained at the high 
end of the singlet mass range shown in Fig. [^ when an¬ 
nihilations into muon pairs are kinematically accessible. 
In that case, F^)^" (m = 2ms) is given to a good approx¬ 
imation by the decay width into muon pairs. 


r{h* gg) = ^—^{2ms) 

Stt 


Inserting this into Eq. (27 ) yields 


avre\ cs (1.5 X 10"“ GeV"®) A® 


m, 

1 - - 

mi: 


1 - 


-1 3/2 


-1 3/2 


(28) 


(29) 
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Even ignoring the kinematic suppression from the muon 
mass and the upper bound on |Ap| from LHC Higgs data, 
it is not possible to obtain a large enough dark mat¬ 
ter annihilation cross section while satisfying the con¬ 
straint on Xp from perturbative unitarity; one would need 
|Ap| ~ Stt, in violation of Eq. (12). Imposing the con¬ 
straint on |Ap| from LHC searches for invisible Higgs de¬ 
cays or the measurement of the Higgs signal strengths 
yields an annihilation cross section more than six or¬ 
ders of magnitude too small to yield the correct thermal 
freeze-out dark matter relic abundance. 

We also note that, even if the appropriate annihilation 
cross section for thermal freeze-out could be obtained, 
this low-mass scenario with SS —)■ is strongly ex¬ 

cluded by cosmic microwave background constraints aris¬ 
ing from the electromagnetic energy injection at the time 
of recombination [?S] (for an update with Planck data see 

Ref. [SO]). 


3. Relic density: thermal freeze-in 

The singlet scalar dark matter extension of the SM 
can easily account for the observed dark matter self in¬ 
teractions if the singlet mass is below about 0.1 GeV, as 
shown in Fig. The model fails when it attempts to 
account for the dark matter relic abundance via thermal 
freeze-out, because the singlet annihilation cross section 
into SM particles cannot be made large enough. How¬ 
ever, there is an alternative perfectly viable way to obtain 
the correct dark matter relic abundance through thermal 
freeze-in nisi]. In this scenario, the cross section for 
two dark matter particles to annihilate into two SM par¬ 
ticles, or vice versa, is extremely small, so that the dark 
matter never comes into thermal equilibrium with the 
SM plasma in the early universe. Instead, dark matter 
is slowly produced through annihilations of SM particles, 
until the universe becomes too cool for dark matter to be 
produced. The resulting dark matter abundance “freezes 
in” and persists to the present day. Because annihila¬ 
tion of dark matter particles never becomes important, 
the kinematic thresholds shown in Fig. [^ are irrelevant 
in this scenario. 

This scenario was studied for the singlet scalar exten¬ 
sion of the SM in Refs. [HHSj. In our scenario with 
^ 0-1 GeV, the correct dark matter relic abundance 
is obtained for |Ap| ^ (2 — 3) x 10“^^ [ISlQ easily 

consistent with the upper bounds |Ap| < (6.1 —9.3) x 10“^ 
from LHG Higgs signal strength measurements and invis¬ 
ible decay branching ratio limits. Such a small 5'5'$1<I> 
interaction strength is even technically natural: if set to 
zero, Ap cannot be regenerated through loops involving 
other couplings in the model. The only disadvantage is 


^ Note that, for this coupling strength, the ApSS<t>t(J) term con¬ 
tributes about (1 MeV)^ to 11311511 . 


phenomenological: the dark matter becomes impossible 
to detect through non-gravitational interactions with SM 
particles. 


III. SINGLET EXTENSION OF THE TWO 
HIGGS DOUBLET MODEL 

To try to salvage the thermal freeze-out scenario for 
singlet scalar dark matter, we consider extended models 
with two resonances, the discovered SM-like Higgs boson 
h with mass 125 GeV and a second CP-even scalar H, 
which can be lighter or heavier than h. Generic models 
allow both states to couple to SS; as a prototype we 
consider the singlet scalar dark matter extension of the 
two-Higgs-doublet model (2HDM) [32lI42| . 

As mentioned previously, thermal freeze-out in the low- 
mass scenario with ms ^ 0.1 GeV is already strongly ex¬ 
cluded by cosmic microwave background constraints [49] . 
We therefore focus in what follows on resonant scattering 
and annihilation through the second scalar resonance H. 
In the singlet scalar extension of the SM, the resonant¬ 
scattering scenario failed because the cross section for 
dark matter annihilation via the Higgs resonance was 
several orders of magnitude too large. In the singlet 
scalar extension of the 2HDM, we can get around this 
by having the dark matter interact resonantly through 
the H pole while making the coupling of H to SM par¬ 
ticles sufficiently weak, thereby suppressing the thermal 
dark matter annihilation cross section through the H res¬ 
onance. As we will see, this will not work in the usual 
2HDMs with natural flavor conservation, but it can be 
made to work in the Yukawa-aligned 2HDM extended 
with a singlet. 

Of course, the singlet scalar extension of the 2HDM 
can also account for the observed dark matter relic den¬ 
sity through thermal freeze-in when ms (5) 0.1 GeV and 
the couplings of S to the two Higgs doublets are both 
extremely small. 


A. 2HDM with natural flavor conservation 

The scalar potential for the singlet scalar dark matter 
extension of the 2HDM can be written as 

V=^s^ + XsS'^ + Api.5.5$l$i + Ap2^^$|$2 + E2HDM, 

(30) 

where <i)i and $2 are the two Higgs doublets and V 2 HDM 
is the usual scalar potential of the 2HDM. We have im¬ 
posed a second Z 2 symmetry under which $1 —>■ — $ 1 , 
which can be softly broken by a dimension-two term in 
V 2 HDM- This allows us to avoid flavor-changing neu¬ 
tral Higgs interactions through natural flavor conserva¬ 
tion |52| . by forcing the right-handed fermions to couple 
to only one Higgs doublet. The transformation proper¬ 
ties of the right-handed fermions under the second Z 2 
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symmetry determine the “Type” of 2HDM (for a review 
see Ref. [53]). 

To set the notation, we write the vacuum expectation 
values as 



with vf + V 2 = v'^ (246 GeV)^ and tan/3 = V 2 1vi. In 

terms of the real neutral components of the two doublets 
4)9 = 0°’”/-\/2, the physical CP-even neutral scalars h 
(identified as the 125 GeV Higgs) and H are defined as 

h = —sma(Pi -I-cos a 02 , 

= cosq; 0°’’'-I-sina^^’*^, (32) 


where a is a mixing angle. 

In terms of these parameters, the Feynman rules for the 
SSh and SSH vertices are igssh and igssH, respectively, 
with 


gssh = 2w(Aj,i cos (3 sin a — Xp 2 sin/3 cos a), 
gsSH = 2u(—Api cos /3 cos a — Xp 2 sin /3 sin a). (33) 

Perturbative unitarity of the SS —>■ 4>i4)i and SS 
$ 24)2 scattering amplitudes separately constrair 0 


|Api|<47r, |Ap2|<47r, (34) 


so Igssffl can be as large as 8ttv. 

In the 2HDM plus a singlet scalar, the low-energy ef¬ 
fective SS —>■ SS coupling for ms mh, mn is given by 


Aeff = As 


9ssh 

8ml 


9ssh 

8m'jj 


(35) 


The constraint ms ^ 0.1 GeV carries through from the 
analysis in the previous section, so we do not pursue this 
low-mass scenario further. 

To avoid experimental constraints from h coupling 
measurements, we tune sin(/3 — a) ~ 1. This means that 


sin Of cos 0, cosQ;~sin/3. (36) 

For sin(0 — a) ~ 1, the H coupling to vector boson pairs 
is heavily suppressed, being proportional to cos(/3 — a). 
This allows H to evade searches for invisibly-decaying 
Higgs bosons via H production in vector boson fusion 
or associated VH production. We focus on the Type- 
I 2HDM, since in that model all the fermion Yukawa 
couplings can be suppressed by making cot j3 small. The 
H couplings to fermions in the Type-I 2HDM are equal 


® A proper coupled-channel analysis would lead to a tighter con¬ 
straint. 


to the corresponding couplings of the SM Higgs boson 
multiplied by a common scaling factor 


sin a 

K 7 = ^~ - cot / 
sm/3 


V 


(37) 


where we used Eq. (36 1 to obtain the coupling in the 
sin(/3 — a) —>■ 1 limit. Perturbativity of the top quark 
Yukawa coupling requires cot/3 <10/3 [51] . 

For ms ~ mij/2, the SS —)■ SS scattering is domi¬ 
nated by the s-channel H exchange diagram, yielding a 
cross section per unit dark matter mass 


As in the case of the singlet scalar extension of the SM, 
the singlet mass must be tuned to be extremely close to 
mn /2 for this scattering cross section to account for the 
halo displacement observed in Ref. |5]. 

The dark matter annihilation cross section to SM par¬ 
ticles, which controls the relic density, is similarly given 
for annihilation through the H resonance by 




0^2 

^ySSH^ tot_ 

EcmiE^m ~ 


(39) 


These expressions take an analogous form to those in 
Eqs. 0 and ( [m] ) for the singlet scalar extension of the 
SM. In that model, the problem was that the annihila¬ 
tion cross section was much too large once the dark mat¬ 
ter self-interaction cross section was fixed to the newly- 
observed value. We therefore attempt to suppress the 
width of H into SM particles while keeping gssH largej^ 
We take mn > nih, so that ms > mhl‘2. and there is 
no need to tune gssh to be small, since invisible decays 
h —>■ SS will be kinematically forbidden. To suppress 
H —>■ WW, ZZ and keep the h couplings SM-like, we re¬ 
quire sin(/3 — a) —>■ 1. To suppress H decays to fermions, 
we take cot /3 1 in the Type-I 2HDM. Then the de¬ 

cay width of H to SM particles, for mn < 2mt, can be 
approximated as 


(42) 

where BR(/ism ^ VV) = BR(/ism ^ WW) + 
BR(/isM —t ZZ) ~ 0.24 |35|. In particular, oc 

(k^)^ ~ cot^ /3. 


Note that, if mu < the h SS invisible decay becomes 
kinematically accessible and must be avoided by tuning gssh to 
be small. Together with Eq. | |36[ |, this implies that 

Api cot/3 ~ —Ap 2 tan/3, (40) 

and hence 

— 2nApi cot^ ~-2uAp2 tan^. (41) 
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Comparing Eq. (14) with Eqs. (39) and (42), we see 
that we need to severely suppress the parameter combi¬ 
nation {nf)^/gggfj. However, in this model the SSH 
coupling is 


9SSH - 2vsmPcos/3{-Xpi + Xp2), (43) 

where we have used sin(/3 — a) ~ 1. In particular, 
dSSH oc. cos/3, so that the ratio / Ossh 

all suppressed at cot/3 1. Indeed, once the upper 
bounds |Api| < 47r and |Ap 2 | < 47r from perturbative uni- 
tarity are imposed, the situation is no better than in the 
singlet scalar extension of the SM discussed in Sec. m 


B. Yukawa-aligned 2HDM 


In the preceding section our main problem was that 
suppressing the H couplings to fermions simultane¬ 
ously suppressed the SSH coupling, so that the fac¬ 
tor /gggfj in Eq. (39) remained no smaller than 

in the singlet scalar extension of the SM. This could 
be avoided if the scalar potential contained a term 
(S'S '$|$2 -I- h.c.). Such a term breaks the Z 2 symmetry 
imposed to avoid flavor-changing neutral Higgs interac¬ 
tions; however, there is no reason to assume a priori that 
flavor conservation originates from the structure of the 
scalar potential. 

The Yukawa-aligned 2HDM allows both Higgs 
doublets to couple to all fermions, and avoids flavor¬ 
changing neutral Higgs interactions by positing that the 
(unknown) flavor physics causes the two Yukawa matrices 
for the up-type quarks to be proportional to each other, 
and similarly for the down-type quarks and the charged 
leptons. Then one can add a third term to the scalar 
potential in Eq. (30) coupling the singlet scalar to the 
doublets. 


V D (Ap3S'5'$l$2 -b h.c.), (44) 


where for simplici^ 
ing Ap 3 to be realn 


we will avoid CP violation by choos- 
The couplings in Eq. (33) become 


gssh = 2u[Api cos /3 sin a — Xp 2 sin /3 cos a 
—Ap3(cos/3 cosa — sin /3 sin a)], 
gsSH = 2u[—Api cos j3 cos a — Xp 2 sin /3 sin a 

—Ap3(cos/3sinQ;-b sin/3cosa)]. (45) 


Tuning sin(/3 —a) ~ 1 to make the h couplings SM-like 
reduces these to 


gssh ^ 2v[-Xpi cos^ P - Xp 2 sin^ l3 - 2Xp3 sin/3 cos/3], 
gssH ^ 2w[(-Api -b Ap2)sin/3cos/3 

—Ap3(sin^ P — cos^ /3)]. (46) 


Now a suppression of the factor in Eq. (39) need 

not be counteracted by a suppression of gggH- We can 
choose the relative Yukawa couplings of the two Higgs 
doublets so that the H couplings to fermions are all 
suppressed while at the same time keeping gssH un¬ 
suppressed. We need a fermionic coupling suppression 
kJ ~ 10“"^ in order to obtain an annihilation cross sec¬ 
tion in the right range for thermal freeze-out. For a full 
treatment one should also take into account the effect 
of thermal averaging of the annihilation cross section at 
the time of freeze-out, which may be important due to 
the severe tuning of the singlet mass relative to the H 
resonance. We leave this for future work. 


We finally remark that, if mn > ruh, this scenario of¬ 
fers the possibility of direct dark matter detection due 
to scattering through exchange of h. The SSh coupling 
need not be small if /i —>■ SS is kinematically forbidden. 
However, the couplings of H to fermions (and hence to 
the nucleons) must be severely suppressed. Thus the di¬ 
rect detection cross section in this scenario is controlled 
by different model parameters than the annihilation cross 
section, so the two are not necessarily correlated. There¬ 
fore a direct detection signal in this model would not 
necessarily probe the coupling responsible for the relic 
density. 


IV. CONCLUSIONS 


In this paper we use the apparent observation of dark 
matter self-interactions in the cluster Abell 3827 to con¬ 
strain models of singlet scalar dark matter. We showed 
that the most natural scenario in this class of models 
is very light dark matter, below about 0.1 GeV, whose 
relic abundance is set by freeze-in through extraordinar¬ 
ily weak interactions with the Higgs doublet (s) through 
a quartic coupling of order 10“^^. 

Thermal freeze-out can be salvaged by extending the 
Higgs sector with another doublet, so long as flavor con¬ 
servation is achieved through Yukawa alignment. The 
solution is, however, quite fine tuned: the singlet mass 
must be tuned very close (of order a part per mille) to 
half the mass of the second CP-even neutral Higgs boson. 
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